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Оптимизация метода Питасси 
вычисления свертки 


Предложенный метод расширяет диапазон используемых разрядностей циклической свертки за счет 
применения эффективного метода вычисления циклической свертки разрядностью 2К, где К — нечетное. 
Показано, что для вычисления свертки такой разрядности достаточно вычислить только две свертки 
половиной (от начальной) разрядности, при большем количестве пред- и поствычислений в виде 
циклических сдвигов по сравнению с методом Питасси. Представлены в общем виде формулы 
вычисления циклической свертки. Приведена реализация операции многоразрядного умножения на 
основе циклической свертки. В виде таблицы приведены оценки сложности вычисления свертки 
большой разрядности вида М=К.2", п>1 для К = 3,5,7,9. 


Введение 


В последнее время большее внимание уделяется вопросам цифровой обработки 
сигналов и методам оптимизации вычислений, связанных с этой обработкой [1-6]. 
В статье особое внимание уделяется методам вычисления циклической свертки, которая 
широко применяется в построении цифровых фильтров, асимметричной криптографии и т.д. 
Как известно, выбор метода зависит от его области эффективного использования. Цик- 
лическая свертка применяется при реализации быстрой операции умножения. Операция 
умножения занимает большую часть вычислительного времени в операциях асим- 
метричной криптографии, таких, как генерация и распределение секретного ключа, 
шифрование и дешифрование информации, проверка и наложение цифровой подпи- 
сиит.д. 


Целью данной статьи является оптимизация вычисления свертки. Основное 
внимание уделяется уменьшению сложности алгоритма за счет уменьшения числа 
вычислений сверток меньшей разрядности. В данной работе приводится применение 
циклической свертки для реализации операции умножения. Как известно, эффектив- 
ность операции умножения больших чисел зависит от метода реализации, а точнее, 
от количества операций умножения однократной точности, поэтому основной целью 
статьи является уменьшение общего числа операций умножений. 

Ниже описан метод, при котором для вычисления свертки разрядностью 2К, 
где К - нечетное, достаточно вычислить 2 свертки разрядностью К. Приведем вна- 
чале эффективный алгоритм вычисления свертки разрядностью 2К, где К -— четное, 
предложенный Питасси, который требует вычисления 3 сверток меньшей разряд- 
ности на всех итерациях. 
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Операторы и условные обозначения 


№ У 
Последовательности разрядности № могут быть представ- зт_| М =: 2 
` Е . > ый . 
лены в виде векторов: в 
Хм Ум- 


Циклическую свертку двух последовательностей х и у обозначим оператором ©: 
М-—1 Е 
Фу РГЕЩЬАЬЕ на ‚ Е=0, М-1, 
р=0 


что равносильно записи: 


ж% | У Я № ое а 
х|я № № В а 
ХУ У № У о о 
х | Ж% Я У а а М ВЫ 

7” 3 Л В п п г В 


Далее будет использоваться именно это представление. 


Рассматривается циклическая (или еще известная как арифметическая) свертка, 
поэтому слово «циклическая» в дальнейшем опускается. 


Операторы Е,О,Г,И (Е-Еуеп (четный), О-ОЧА (нечетный), [.-Гомег (нижний), 
О-Оррег (верхний)) определим следующим образом: 


(Ех) =», (ОХ), = хи ь (1%), = х, ‚ (ИХ); = Хуурь› К =0,№/2-—1. 


! 


Циклический сдвиг влево элементов обозначим через 2 = (у ) : 


2. = Уцчутом › 20, М-1. 


т 


Циклический сдвиг вправо элементов обозначим через 2 = (у ) : 


2, = Ум шоа м > [4 = 0,М-1 я 


Проиллюстрируем вид операторов для случая №М=8: 


№ Хх Хо Хх У Ув 
Хх Хх Хх Хх ' " 
| С бы 3 = 1 | в: | (Е) = У (Е > У 
х. х х> % У У> 
% Х' Хз Х' У У. 


уп 
Циклический сдвиг влево элементов и раз обозначим через 2 = (© ) } . 


п\п 
Циклический сдвиг вправо элементов п раз обозначим через 2 = (© ) } р 
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Метод парисекции 


Рассмотрим метод, предложенный Питасси [3] и обобщенный Девисом [6]. 

Лемма 1 [3]. Входящие и исходящие последовательности свертки Г =х®у могут 
быть расщеплены на равные группы в соответствии с четностью элементов последователь- 
ностей: 


Ет = Ех ® Еу + Ох ® Оу, 
От = Ех® Оу+Ох®(ЕУ). 


ЗЕ М 3 р м о в 
м № м м Хх Хх № м У | м № №№ № м 
о ее Е еды Па Мо 6 


т Е а 
а” Же О Е 9% ПЕ; Обе Эа О ей о» 96 
о о Ве О. В. В. Е о с 
ЕТ Ог 
Рисунок 1 — Вычисление свертки Рисунок 2 — Вычисление свертки 
для № =8 стандартным методом для М =8 методом парисекции 


Из выражений: 
а =(Ех + Ох) ® (Еу+ Оу), 


Ь =(Ех-0)®(Е’- ОУ), Е" =У2(а+5), 
с-0т® (Е) -Е5 Ог=12(а-Б)+с 


можно получить соотношения (1). С более детальным выводом формул (2) можно 
ознакомиться в [3]. 

Перейдем к вычислению свертки разрядностью М =2К, где К - нечетное. 
Рассмотрим пример свертки разрядностью М =6=2.3: 


(2) 


У № № № № № М Ум м №№ № М 
Ям № № м № Х Еу У, | Хх Хх Жим Хх ХХ: 
Ул М № № № И В И Е 
У № %Х ХМ № № № Ям № №№ м ХХ 
У | м № № м м ХХ Ом м мих м 
У № № № № м № У № № №1№м № № 
п В Ц ББ п р вап В п Б 
Е’ О’ 
Рисунок 3 — Вычисление свертки Рисунок 4 — Вычисление свертки 
для М№ = 6 стандартным методом для М№ = 6 методом парисекции 
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Воспользуемся тем свойством свертки нечетной разрядности, что циклический 
сдвиг нечетное число раз в одном направлении дает результат четного сдвига в дру- 


гом направлении. 


а = | В+ (0%) |® 


Ы + (05) | 


„- 


= | Е - (0%) |®| 7-() | (3) 
Е = 29 +2)= Ех ® Еу+ (0%) ®(05) = Е® +0х®0у, (4) 
о ем | Ее ® (0) +©) ФУ = ® 0+0 ®(Е). (5) 
Лемма 2. Имеют место соотношения: 
(0х) ® (07 =0®0у, 
Е; ® (0) = (Е ®0), (6) 
ох®(Е>) = (0х) 85. (7) 
Рассмотрим последовательно каждое соотношение на примере М =6. 
Начнем с (Ох) © (ОУ) =Ох®оду. 
Ух Хх Хз Ух Хы ХХ 
Оу ух хх х; (0У) У 1х; № Я 
У; Хх; Хз Хх У! Хх Х5 Хз 
20 2. 24 20 2. 24 
Ей Е? 
РО ЕЕ = (0%) ® (0) 
Теперь рассмотрим Ех © (0) = (Ех © 0) з 
Я | № 4% 4% | У | № и № 
ОУ ух жЖ Хх, (0у) У; | № № х4 
У |4 № № Я | %4 х № 
2 С С: Язь 255 21 
о (02) 
О; = Е ®Оу (©=) - = ® (05) 
Аналогично можно показать, что Ох © (Е). = (Ох © ВУ). р 
Из (6) следуст, что: [2 ® (©) а. (Е ® сз) - Х®Оу. 
В заключении рассмотрим Ох ® (Е). = (©) © Ру. 
Ух ХХ; ФХ У | м Х5 Хз У | 23 х! № 
Еу Ум а (Е) Зе | х № ЕУ У, | х; №3, ый 
Ух; № Хх У. |х. х, х, У4 | Хх. Хз 
2 2; 25 25 2 2 25 р в 
О (0=)' (07) 
ОЕ-ОЕ® Е (02). = ®(Е. (0=)' = (0х) ®Е 
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Аналогично можно показать, что верно соотношение (Ох) © Еу=Оу® (Ех) . 
С учетом формул (6), (7) соотношение (5) примет вид: 


[а 


(©) 25) = (0х) ®(ЕУ) = 0х (5) (8) 


С учетом того, что два циклических сдвига вправо можно заменить одним 
сдвигом влево, соотношение в формуле (8) примет вид: 


Ох® (5) - 0х8 (5). 


Продемонстрируем более подробно, в чем состоит предложенный метод. 


м м №1№ № М Хм м №№ х № 
Е У | х № Жжах № № Е Ух № Хх № х 
о а Е В Ям № мм № № 
м № №№ № № | Я М %\% м м № 
О ух м Мих м Хх (0У) У № № № м № 
У | № № №10 ХХ № Ям № №1м ^ № 
С Е. с. В Е АЕ. В я 

р 

ЕР О’ ЕР | (ОР) 

Рисунок 5 — Вычисление свертки Рисунок 6 — Вычисление свертки 
методом парисекции предложенным методом 


Из рисунков видно, что результат совпадает. Так, например, при использова- 
нии обоих методов и равно: 
7= Ух + Ум + УХ + Ух, + УЖ т УЖ, 


= Хо + МУ, + ХУ, + ХУ. + Х.У; ХУ, 


5 
п= Уд Ус) шоаб * 
0 


Видно, что в методе парисекции (рис. 5) нужно делать циклический сдвиг толь- 


и 
ко в векторе (Ох) ‚ тогда как в предложенном методе, согласно соотношению (3), 


1 1 
нужно циклически сдвигать в противоположную сторону оба вектора (Ох ) . (ОУ) перед 
вычислением свертки. Из рис. 6 также видно, что для получения корректного резуль- 


т 


тата нужно сдвинуть вектор (© ) } =О!. 


Предложенный метод требует больше пред- и поствычислений в виде цикли- 
ческих сдвигов, но позволяет вычислить свертку за счет только 2 сверток меньшей 
размерности вместо 3 согласно стандартному методу парисекции. 

Достаточно 10 умножений для вычисления свертки разрядностью М№М=5 [2] и 
соответственно 20 умножений для вычисления свертки разрядностью № =10 =2.5. 
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В этом случае формулы (3) - (5) примут вид: 


р в+(©) ® в+(©5) 


== ве-((0) ® в-(©) | 


ЕР=2( +2), 


От= уз(@ -=} | 


Из формул видно, что нужно сдвинуть нечетные элементы на 2 позиции. 


У | № 4 № 4 | № м ХХ № М 

У | ХХ № № М | мо № Ху № М 

Е У. | х4 Хх № ЖЖ | о. 
Ув | № м № № № № № М № Ху 

Уз | № № № ХХ. № | ХХ; Хз М _ № 

У № м № Хз | Хо № № №4 ХХ 

У Л № Хр № | Хо № ХХ  Х4 

Оу ух Хх № Х7 | д М № № № 
У хо о № №4 № № ЖЖ 

У | № № ХХ 2 м | № № №4 №2 № 
Г) й РР | в И #: Ри 

ЕР ОГ 


Рисунок 7 — Вычисление свертки разрядностью М =10 методом парисекции 
Продемонстрируем, что, используя всего две матрицы меньшей разрядности, 


ое о о о 
И а. Е дом № № 
В. ЗФ ие СА (©) ] Е 
т ее Ж № № 
Г а, о Е ба 508 
можно представить свертку разрядностью № =109 =2.5: 
У | № № № № 4 | № м № 
| № № № | № № Хз м № 
Еу Ум № № № ЖЖ х5 3: М 
У | № № № № м № Х7 5 № 
_ У №64 №2 № | м _ №. №7 _ №5 
5 Из ЖА ОХ | № № № № 
Ум ож м м № № № № № 
‚2 | 
(>) Ум м о Х м мм № ж № 
| 
Ям № м № № | № №4 Х> о: 
У | м № № № № Х4 Хх. № 
п Бо К В | 5 п ИЕ В. 
2 
ЕР | (=) 
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Количество сдвигов равно и =| К/2 |, где К - нечетное число и множитель 


разрядности свертки №М=2К. 


Теперь можно выразить в общем виде формулы для вычисления свертки 
разрядностью М№ = 2К ‚ где К - нечетное число. 


а 


[Е й (©) ] ® | х (©) | 


и [= (5) |®| 5-(©)) | 


ЕЕ = 2+2), ОЕ = уз(@-е] | | 


Покажем практическое применение свертки для вычисления операции умножения 
многоразрядных чисел, которая широко используется при реализации алгоритмов асим- 
метричной криптографии. Покажем вначале на примере свертки разрядностью 2К, 


где К`— четное. 


Умножение двух чисел х = в = (х, Уы У У.) разрядностью 4 с ис- 


пользованием свертки можно представить в виде: 


в = (хх, х,, №,0,0,0,0)® (0,0,0,0, ууу, У). 


х |0 0 0 Оу У Ух № 
х |0 0 Оу ил 0 
х, |0 О у ил 0 0 
х |0 у У У 0 
Оу у у р 0 0 
Оу хр 0 0 0 0 у, 
Оу у 0 0 0 0 у; у, 
Оу 0 0 0 Оз у У 
о 16. о а С С 


Рисунок 8 — Вычисление свертки 
стандартным методом при реализации 
операции умножения двух чисел 
разрядностью 4 


Окончательно получаем: 


№ 33: 92. Е 6 
м Эа дих И 520 
^2 23. Ир М 20 
№ УЗ. а 12 
с ЗА. В В 


Рисунок 10 — Умножение двух чисел 
с использование свертки без нулевых 
элементов 
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0 0 0 О у У № 


Оооо У ул 0 
Обузу у 0 0 
Орли 0 0 0 
оо оооо о о 
оо ооооо о 
оо ооооо о 
ооо оооо о 


7 Е Г 7. Г, Г й т 


Рисунок 9 — Вычисление свертки 


стандартным методом при реализации 


операции умножения двух чисел 
разрядностью 4 с учетом нулевых 


элементов 
Уз У> У Ус 
Хх х› х № 
ХУ №, У №У 

ЖУ МУ МУ МУ 
ХУ ХУ №У ХУ 

ХУ ХУ №) ХУ 

” 7% Г 74 7 7) п о 


Рисунок 11 - Умножение двух 


чисел в столбик 
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Из рис. 10 и 11 видно, что использованная таким образом свертка реализует 
стандартный метод умножения в столбик, за исключением того, что результат 
получается в обратном порядке. 


Приведем пошаговое описание алгоритма реализации операции умножения К- 
разрядных чисел с использованием свертки: 

Шаг 1.а; =х,, 1=0,К-Т,а, =0,1=К, М-1. 

Шаг 2. 5 =0,1=0,К 1,5 =у 1, Ё=К,М-Т. 

Шаг 3.7 =а®Ь. 

Шаг 4. \, =г 


ЕР, Ё= О, М-1. 

Вектор м разрядностью № =2К будет содержать результат операции умно- 
жения чисел, представленных в виде вектров х и у разрядностью К. 

Данный алгоритм справедлив и для сверток рязрядностью М =2К, где К — 
нечетное. 


Заметим, что для вычисления операции умножения К-разрядных чисел с 
использованием алгоритма быстрого преобразования Фурье необходимо вычислять 
циклическую свертку также с разрядностью 2К. Это связано с тем, что при умно- 
жении двух К-разрядных чисел, результат умножения получается разрядностью 2К. 


Количество операций умножения при вычислении свертки разрядностью 
№ = 2”, п>2 методом парисекции выражается следующим соотношением: 


Му =5.3=2, (9) 
где М - количество операций умножения, Л! — разрядность свертки. 


Количество операций умножения М при вычислении свертки большой 


разрядности вида М№М=К.2”, п>1, где К - нечетное, приведено ниже в виде 
таблицы. 


Таблица 1 
Количество операций | Разрядность Количество операций 
Значение | умножения М „ для свертки умножения ЛИ при 
Кв (9) вычисления свертки А. вычислении свертки 
разрядности К п>1 разрядности № = К.2",п>1 
3 ь №М=3.2" М =8.3"" 
5 10 А М =20.3"" 
Г 16 др" М =32.3"" 
2 19 №=9.2" М =38.3"" 
Из табл. 1 видно, что можно конструировать свертки размерностью, отличной 


от М=2”", и>2, что является одним из преимуществ предложенного метода. Вто- 
рым преимуществом является то, что в некоторых случаях эффективнее использо- 


вать свертку большей разрядности вида М№М=К.2”, п>1, где К - нечетное, вместо 
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свертки разрядностью М№М=2", и>2. Так, например, для вычисления свертки раз- 
рядности № =18 требуется всего 38 умножений согласно таблице, тогда как для 
вычисления свертки М№=16 методом парисекции необходимо 45 умножений по 
формуле (9). 


Заключение 


В данной работе приведен эффективный метод вычисления свертки разрядностью 
2К, где К - нечетное. В общем виде приведены формулы вычисления циклической свертки 
разрядностью 2К; где К’ нечетное. Приведен алгоритм реализации операции умножения с 
использованием циклической свертки. В виде таблицы приведены оценки сложности 


вычисления свертки большой разрядности вида № = К-2”, п>1, где К = 3,5, 7,9. 
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ОптимзащЯя методу Штасс! обчислення згортки 

Запропонований метод розширюе д!апазон використовуваних розрядностей цикл!чно! згортки за 
рахунок застосування ефективного методу обчислення цикл!чно! згортки розряднастю 2К, де К -— 
непарне. Показано, що для обчислення згортки тако! розрядност! достатньо обчислити Пльки дв! 
згортки половинно! (в1д начально!) розрядность, при б1лышй клькост! перед- та постобчислень у 
вигляд! цикл!чних зсув1в. Представлен! в загальному вигляд! формули обчислення цикл!чно! згортки. 
Наведена реалзацй операцИ багаторозрядного множення на основ! цикл!чно! згортки. У вигляд! 
таблиц! наведен! оцнки складност! обчислення цикл!чно! згортки велико! розрядност! виду М=К.2”, 
п> 1 для К = 3,5,7,9. 


А.М№. Тегезйсепко 

Тре засое$е тео ехеп45 Фе гапое оЁ изе4 теазигетлеп оЁ сусПс сопуойоп$ У ито оЁ еЙеснуе 
сайсщаноп тефо4 оРсусПс сопуо[Нопз ул теазигетлепе 2К еп К 1$ 044. [1$ $ЗПохуп Юг сопуоайоп са1сиайоп 
ул Фа еазигеглепе # 15 епоизВ {0 сасшае ошу 2 сопуоавоп$ Ва тлеазигетепе ул тоге пипабег оЁрге- ап 
роз{-са1сайоп$ ПКе сусПс $015. [4 1$ отуеп 1 оепега] ве сасшайоп Югицаз оЁ сусПс сопуо[айоп. Ё 1$ этуеп фе 
Байате ог та@-Фей ши@рИсайоп мб чзше сусПс сопуоайоп. Тре сотар!ехез оЁ сусПс сопуоайоп 
са1сщавоп уу плеазигетепЕ № = К’.2”", и>1 Юг К = 3,5,7,9. аге &1уеп ш 1аЫе. 


Статья поступила в редакцию 21.07.2008. 
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